Fast population transfer engineering of three-level systems 
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We design, by invariant-based inverse engineering, resonant laser pulses to perform fast popula- 
tion transfers in three-level systems. The efficiency and laser intensities are examined for different 
protocols. The energy cost to improve the fidelity is quantified. The laser intensities can be reduced 
by populating the intermediate state and by multi-mode driving. 
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I. INTRODUCTION 



The laser control of internal state preparation and dy- 
namics is of importance in atomic and molecular physics 
for applications such as metrology, interferometry, quan- 
tum information processing and driving of chemical re- 
actions [H-§. In two- or three-level systems, resonant 
pulses, rapid adiabatic passage (RAP), stimulated Ra- 
man adiabatic passage (STIRAP), and their variants 
have been widely used to perform population transfers 
0-3- Generally, resonant pulses may be fast, but are 
highly sensitive to the deviations of pulse areas and exact 
resonances, whereas the adiabatic passage techniques are 
robust versus variations in experimental parameters but 
slow. To combine the advantages of resonant pulses and 
adiabatic techniques and achieve fast and high-fidelity 
quantum state control, some alternative approaches, like 
composite pulses and optimal control theory (9Hllj , 
have been proposed. Several recent works on "shortcuts 
to adiabaticity" have been also devoted to internal state 
population transfer and control (T^ - [l~9| . The shortcut 
techniques include counter-diabatic control protocols jl2l 
or, equivalently, transitionless quantum driving jl3l4l5l |. 
"fast- forward" scaling [2(J , and inverse engineering based 
on Lewis- Riesenfeld invariants pll . l22j . These methods 
are in fact strongly related, and even potentially equiva- 
lent (lol l2p. However, they provide in general different 
shortcuts jl6l. Il9j . 

In this paper, we apply invariant-based engineering to 
realize fast and robust population transfers in three-level 
systems. This method has been applied in trap expan- 
sions ^mHiHiz!], rotations [H], atom transport [H- 

t mechanical oscillators [lH, or many-body systems 
|34| . In a three- level system as the one depicted in 
Fig. [TJ STIRAP allows to transfer the population adia- 
batically from the initial state |1) to the target state |3). 
To speed up the process, a fast-driving counterdiabatic 
field connecting levels |1) and |3) may be used [l"H [35|. 
In general, though, this implies a weak magnetic dipole 
transition, which limits the ability of the counterdiabatic 
field to shorten the times [HI, [35j . This will be solved by 
invariant-based engineering, which provides alternative 
shortcuts without coupling directly levels |1) and |3). 
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FIG. I: Level scheme of STIRAP for a A level configuration. 
fi p and Q s are the Rabi frequencies for the interactions with 
the pump and Stokes fields respectively, and A p and A3 are 
the detunings from the resonances. 



II. INVARIANT DYNAMICS 

To perform STIRAP, the Hamiltonian within the ro- 
tating wave approximation (RWA) reads 

* / a p (t) \ 
H (t) = o n P (t) 2A P n s (t) , (i) 
2 \ n,(t) 2A 3 ) 

where, as shown in Fig. [I] the coupling strengths be- 
tween the states are determined by fl s (t) and £l p (t), 
which describe the interactions with the pump and Stokes 
fields, and the detunings from resonance arc defined as 
A p = (E 2 - EJ/H - w p , A s = (E 2 - E 3 )/h - u a , and 
A3 = A s — A p , where uj s and uj p are the laser (angular) 
frequencies, and the Ej, j = 1,2,3, the bare-basis-state 
energies. 

We consider the so called "one-photon resonance" case, 
Ap = A3 = to simplify the Hamiltonian as 

h ( n p (t) \ 

H(t) = - fip(t) Q s {t) . (2) 

1 \ n s (t) / 

The corresponding instantaneous eigenstatcs \n), with 
eigenvalues Eq — and E± = ±M7/2 with SI = 



2 



'fijj + n 2 s , are 



and 



KM) 
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sin 

|n±(t)) = -j= I ±1 | , (3) 
V 2 \ cos6> 



where tan# = Cl p /Cl s . When the adiabatic condition, 
\9\ <C is satisfied, perfect population transfer from 
the ground state |1) to the excited state |3) can be 
achieved adiabatically along the dark state |no), using 
the counterintuitive pulse order (Stokes before pump). 
In what follows we shall use a dynamical invariant to 
speed up the population transfer. 

To construct the dynamical invariant, the Hamiltonian 
in Eq. (|2|) can be rewritten as (3(| 



(4) 



where Ki , Ko , and K% are angular- momentum operators 
for spin 1 |36j, 



1 
K x = | 1 | , K 2 






1 I ,K a 
1 




that satisfy the commutation relations 

[Ki,K 2 ] = ik 3 , [K 2 ,K 3 ] = iKx, [K 3 , K x ] = iK 2 . (5) 

The Hamiltonian ([2]) possesses SU(2) dynamical sym- 
metry, and an invariant I(t), satisfying the condition 
dl/dt ee dICtVdt + (l/ih)[I(t),H(t)] = 0, can be con- 
structed as [II, H3 



I(t) = -fi (cos7sin/?A'i + cos7COs/3A' 2 + sin 7^3), 

u I cos 7 sin (3 — isin.7 
= -Oo cos 7 sin /3 cos 7 cos f3 | ,(6) 
^ \ i sin 7 cos 7 cos f3 

where fio is an arbitrary constant with units of frequency 
to keep I(t) with dimensions of energy, and the time- 
dependent auxiliary parameters 7 and (3 satisfy the equa- 
tions 



1 



(fi p cos P — Q s sin 



(7) 



/3 = -tan7(Q s cos ) 8 + fipSin/3). (8) 

The eigenstates of the invariant /(i), satisfying 
I(t)\<t>n(t)) = X n \4>n(t)), (we use the labels n = 0,±) 
are 



\Mt)) 



cos 7 cos j3 

—i sin 7 
— cos 7 sin /3 



(9) 



sin 7 cos /3 ± i sin / 
\4>±{t)) = -j= I «cos 7 



(10) 



■ sin 7 sin /3 ± z cos /3 
which correspond to the eigenvalues Aq = and A± 



the 



±1. According to Lewis- Riesenfeld theory 
solution of the Schrddingcr equation, ihdt 1 ^ 
is a superposition of orthonormal "dynamical modes" , 
= En C„e lQ ™|0„(i)) [13, where each C„ is a time- 
independent amplitude and a„ is the Lewis-Riesenfeld 
phase, 

«n(*) = ~ j Q {4>n{t')\lh^ - H{t')\<j> n {t'))dt' . (11) 



In our case ao = 0, whereas 



a± = =F 



/? sin 7 -I — (fJ p sin (3 + fl s cos /?) cos 7 



The dot represents here a time derivative. 



III. INVERSE ENGINEERING AND FAST 
POPULATION TRANSFER 

Beginning with the Hamiltonian in Eq. $2$, we shall 
apply invariant-based inverse engineering to design fl s 
and Q p . We firstly assume that Q s and f2 p are unknown 
functions to be determined, from (JT]) and J5J), 



Q s = 2 ($ cot 7 cos /3 — 7 sin /3), 
f2 p = 2 (/? cot 7 sin (3 + 7 cos /?). 



(12) 
(13) 



Once the appropriate boundary conditions for 7 and (3 
are fixed, we are ready to choose some ansatz, for ex- 
ample, a polynomial or some other function with enough 
free parameters, to construct fl s and Q p . 

Our Hamiltonian H(t), Eq. ([2]), should drive the initial 
state |1) to the target state |3), up to a phase factor, 
along the invariant eigenstate \<fio(t)) in a given time tj. 
We therefore write down the boundary conditions for 7 
and /3, based on Eq. , 



7 (0) = 0, 7 (*/) = 0, 
0(0) =0, £(*/)= tt/2. 



(14) 
(15) 



In general, H{t) does not commute with the invariant 
lit), which means they do not have common eigen- 
states. To achieve fast adiabatic- like passage (i.e., not 
really adiabatic all along but leading to the same fi- 
nal result), one may impose boundary conditions to sat- 
isfy [i?(0),I(0)] = and [H(t f ),I(t f )} = 0, which give 
n p (0) = and fl s (t f ) = 0. Using Eqs. ([T2" l) -([l5 ]) this 
implies the additional boundary conditions 



7(0) =0, 7 (*/) = 0. 



(16) 
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FIG. 2: (Color online) (a) Time evolution of Rabi frequencies, 
fi p (solid red) and Q s (dashed blue) for ~/(t) = e and /3(t) = 
Trt/2tf. (b) Time evolution of the corresponding populations 
of levels 1 (solid red), 2 (dashed blue), and 3 (dotted black). 
Parameters: tf — e = 0.2. 



The set of conditions in Eqs. (fTlT) - p^|) guarantee fast 
adiabatic-likc population transfer. Now we are ready to 
apply inverse engineering and design different protocols. 

Protocol 1- hi the first example, we set the boundary 
conditions for 7 and (3 as follows: 

7(0) = e, 7(0) = 0, j(t f )=e, j(t f ) = 0, (17) 
0(0) = 0, £(*,)= tt/2. (18) 

Note that we have introduced a small value e for 7, as 
an exact zero value implies infinite Rabi-frequencies ac- 
cording to Eqs. (fl"2j) and (fl"3"j) . With these boundary 
conditions, we can simply choose 

7 (t) = e, p(t)=Trt/2t f , (19) 

which provide 

n s (t) = {n/tf)cotecos(nt/2t f ), (20) 



&p(t) — (tt / 1 f) cot e sin (wt / 2t f) . 



(21) 



Fig. [2] shows the time evolution of Rabi frequencies 
and corresponding population transfer for ^(t) with ini- 
tial and final states \<f>o(0)) and \<fio(tf)). We take 
I — 3) = (0,0, —1) T as the target state, which corre- 
sponds to \4>o(tf)) for the ideal conditions j(tf) = and 
0{tf) = (Note that for e ^ the initial state is not 
exactly |1). In Protocol 3, below, we shall examine the 
case 1^(0)) = |1).) The final fidelity with the target state 
is 



FIG. 3: (Color online) (a) Time evolution of Rabi frequencies, 
fi p (solid red) and Sl a (dashed blue), for j(t) = X^=o a J^'' 



and /3(t) = Y2j=o^i^ ' with the boundary conditions (|27|> - 
(|29p . (b) Time evolution of the corresponding populations 
of levels 1 (solid red), 2 (dashed blue), and 3 (dotted black). 
Parameters: 5 — 7r/4, tf = 4/is, e = 0.2. 



From Eqs. l[20] ) -([2"2"] ) . we find 



dfl s dflp ir cos (Trt/2tf) 1 

~ I ~ 2 ~ 2 > ' ) 
t f sm e e 



de de 



and 



OF 
Ik 



(24) 



respectively. In other words, the fidelity varies smoothly 
with e, whereas the Rabi frequencies decrease dramat- 
ically when increasing e. This provides the possibility 
to achieve a desired fidelity with relatively small Rabi 
frequencies. 

Improving the fidelity or shortening 1 1 implies increas- 
ing the Rabi frequencies. Note the behavior of the time- 
averaged frequency, 



7r cot e 



L f Jo v */ 
and the "energy cost" (time-averaged energy), 



E/n= / (cii + n;)dt 

Jo 



tf 



(25) 



(26) 



F 



3|*(*/)) 



cose. 



(22) 



Protocol 2- We design now a different protocol, in 
which the intermediate state |2) may be populated, and 
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FIG. 4: (Color online) Time-averaged frequency (a) and en- 
ergy (b) in the second protocol as a function of <5 for different 
values of e, where e = 0.2 (solid red), e = 0.02 (dashed blue), 
and e = 0.002 (dotted black). Other parameters are the same 
as in Fig. [3] 

both pump and Stokes pulses vanish at t = and t = tj. 
Thus, we set the following boundary conditions: 

7 (0) =e, 7 (0) = 0, j(t f )=e, j(t f ) = 0, (27) 
0(0) =0, P{t f ) =tt/2, (28) 
j(t f /2) = 5, /3(0) = 0, $(t f ) = 0. (29) 

The boundary conditions in Eqs. (f27|) and ([28]) are the 
same as before, but we add now Eq. (|29|): since the 
population of the intermediate state |2) is given by P2 = 
sin 2 7, the condition j(tf/2) — S sets its maximal value 
at t = tf/2, whereas (3(0) = and $(tf) = guarantee 
that il s (0) = and Cl p (t f ) = 0. 

By assuming a polynomial ansatz, 7(f) = X)j=o 

and j3(t) = ^2j=objP , ^° interpolate at intermediate 
times, we can solve the coefficients in terms of the bound- 
ary conditions. Once j(t) and j3(t) are fixed, we may 
calculate the time evolution of pulses and populations, 
see e.g. Fig. [3j where 6 = ir/4 is chosen as an example. 
Fig. |3] shows that the intermediate level |2) is popu- 
lated, and the population is 1/2 at t = tf/2, because 
j (tf/2) = 7r/4. The two examples are compared in Figs. 
[2] and |3l the laser pulse intensity is smaller when the 
intermediate state |2) is allowed to be populated. Note 
that while sharing the SU(2) dynamical symmetry with 
two-level systems [l6|, the three- level system cannot be 
reduced to a two-level system. 

We also calculate the time-averaged Rabi frequency 
and energy in Fig. @] Fig. @] demonstrates that the time- 
averaged frequency and energy increase for a smaller e as 
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FIG. 5: (Color online) Fidelity versus e for the initial state |1): 
solid red and dotted black lines correspond to the examples 1 
and 2 in Figs. [2] and [3] The fidelity F = cose (dashed blue) 
for the initial state |(/>o(0)) is also shown. 



before. When e = 0.002 in Fig. HI the fidelity F = cos e 
is equal to .999998, which satisfies the criterion for a 
fault-tolerant quantum computer [B[ . They also decrease 
significantly by populating level |2), though the behaviors 
of frequency and energy are not the same. Remarkably, 
Fig. |U (a) shows that the time-averaged frequency for 
each e can be minimized. For the smallest e this happens 
when 6 approaches tx/2. In this case, the intermediate 
state is fully populated. The time-averaged energy is 
even flatter for central values of 5, see Fig. |U(b). When 
the intermediate state is not populated at all, that is, 
5 = or 5 = 7T, both time- averaged energy and frequency 
increase dramatically. 

In general we may combine the invariant-based method 
with optimal control to optimize t he p rotocols according 
to different physical criteria [25|, [3(J [3l|, for example, 
(time-averaged) frequency minimization or energy mini- 
mization. The time-optimal problem with bounded en- 
ergy, and the minimum energy cost problem for fixed 
time have been solved for the three-level system || [l(| ■ 

Protocol 3- Our last protocol may be considered as a 
variant of the first one, with the same pulses but a dif- 
ferent initial state. An important difference with respect 
to the previous protocols is that it is based on multi- 
mode driving rather than on a single-mode driving. This 
means that the time dependent wave function ^(t)) will 
include contributions from the three eigenvectors of the 
invariant. 

So far we have assumed that the initial state depends 
on e through the dependence of |0o(O)) on e. Let us 
instead use the bare state |1) as initial state but keep 
the designed interactions e-dependent as before. Fig. [5] 
shows the fidelity ( — 3|^(t/)) as a function of e when the 
initial state is |1). The fidelity for the e-dependent initial 
state |</>o(0)) is also shown for comparison. 

Interestingly, the fidelity oscillates with respect to e for 
the pump and Stokes pulses described by Eqs. (f20|) and 
(f2~T1) . To analyze this in more detail, we first calculate 
the final state V(t f ) = J2 n C«e M " \<f> n {tf)) , where C n = 
(0 n (O)|l). With the eigenvectors \(f> n (t)) at t — and tf 
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FIG. 6: (Color online) Time evolution of the populations of 
levels 1 (solid red), 2 (dashed blue), and 3 (dotted black), 
where the pump and Stokes pulses are described by Eqs. (|20[) 
and (1211) . and the initial state is |1). Parameters: tf = 4^is, 
e = 0.2527. 



we have 



F= (-3|*(t/)) =e lao cos 2 e+-(e M + + e m -)sin 2 e. 



(30) 
i are 

(31) 



In the first protocol, the Lewis-Riesenfcld phases a n are 

7T 

a = 0, a± = Ttt~- — > 
2 sine 

which finally gives 

F - 1 - Sin2e { 1 - COS (^n7)}- 
When the condition 

(sine)" 1 = 4JV, (N = 1, 2, 3...) 



(32) 



(33) 



is satisfied, the fidelity becomes 1. By solving Eq. 
we get e = 0.2527 for 2V = 1, e = 0.1253 for N = 2, etc... 
In particular, the rightmost maximum at e = 0.2527 com- 
bined with the initial state |1) provides stable, perfect 
population transfer, as shown in Fig. [SJ with less intensi- 
ties than the ones required in the first protocol for a good 
fidelity, since the value of e is relatively large now. Com- 
pare the values H = 27r x 0.48 MHz and E/h = 2tt x 5.89 
MHz for e = 0.2527 in Protocol 3 (with fidelity F = 1) 
by using Eqs. fl25) and (JUJ), with = 2tt x 0.62 MHz 
and E/h = 2n x 9.56 MHz for e = 0.2 (corresponding 
to F = 0.9682) in the first protocol. To achieve higher 
fidelity in the first protocol, for example, F = 0.9998, the 
time-averaged frequency and energy cost have to be in- 
creased up to Q = 2?r x 6.25 MHz and E/h = 2tt x 981.49 
MHz by choosing e = 0.02. In summary, Protocol 3 based 
on multi-mode driving provides an alternative shortcut 
to implement a stable, perfect population transfer with 
a low energy cost. 



IV. DISCUSSIONS AND CONCLUSIONS 

In this paper, we have developed the invariant-based 
inverse engineering method to achieve fast population 



transfers in a three-level system. Two different single- 
mode protocols are applied first in which the fidelity is 
linked to the laser intensity. Shortening the time also im- 
plies an energy cost. Interestingly, to achieve the same 
fidelity, less intensity is required when the intermediate 
level 1 2) is populated. A variant of the first protocol in 
which the initial state is simply the bare state |1) and 
the dynamics is driven by a multi-mode wave-function 
provides a stable and less costly shortcut. Further explo- 
ration of the multi-mode approach in this and other sys- 
tems is left for a separate study. The population of the 
intermediate level is usually problematic when its time 
decay scale is smaller than the process time. While this 
may be a serious drawback for an adiabatic slow process, 
it need not be for a fast shortcut. Protocols that popu- 
late level 2 may thus be considered as useful alternatives 
for sufficiently short process times. 

As we stated in the Introduction, different techniques 
to find shortcuts to adiabaticity are strongly related, or 
even equivalent. The invariant-based inverse method pre- 
sented here may be compared to the optimal control ap- 
proach used in Q . In the optimal control method @ , the 
system of control differential equations are the same as 
Eqs. ([7]) and © in the invariant method. The ultimate 
reason is that these equations are in fact equivalent to 
the Schrodinger equation for a given wave-function pa- 
rameterization. The invariant dynamics provide thus a 
complementary understanding of the optimal control ap- 
proach, whereas optimal-control techniques also help to 
optimize the results given by the invariant-based inverse 
engineering. 

Finally, the present results -within the on-resonance 
conditions- are applicable to quantum state transfer with 
three qubits [3^] , adiabatic splitting or transport of atoms 
in a double well, and a triple well [40(. In a more general 
case, the Hamiltonian |T]) (A p ^ and A3 7^ 0) does 
not possess SU(2) symmetry, so that the invariant I(t) 
should be constructed in terms of the eight Gcll-Mann 
matrices for the SU(3) group [4l|. The invariant- based 
inverse engineering for the systems with Gell-Mann dy- 
namic symmetry will be discussed elsewhere. 
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